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ABSTRACT — Russula sichuanensis and its ectomycorrhizae were collected from monoculture 
forests of Pinus wallichiana. This fungus is a new record for Pakistan and the ectomycorrhizae 
are described for the first time, based on morphology and nrDNA ITS sequences. Russula 
sichuanensis is a morphologically variable species especially with respect to the colour 
and shape of its pileus. The taxon is characterized by small to medium-sized stipitate 
basidiomata, reticulate globose to subglobose basidiospores, and pluriseptate pileocystidia. 
The ectomycorrhizae are distinguished by dichotomously branched orange-brown to grayish 
brown tips, its pseudoparenchymatous acystidiate mantle with irregularly shaped and often 
interlocking cells, and curved and frequently septate emanating hyphae. Pinus wallichiana is 
a new host for this species. 
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Introduction 

The genus Russula Pers. is among the most diverse genera of Agaricomycetes. 
The genus includes approximately 750 species (Kirk et al. 2008), many of which 
are notorious for theirintergrading colour variations and obscure morphological 
and anatomical discontinuities (Miller & Buyck 2002). Although thirty-three 
Russula species have been recorded from Pakistan (Ahmad et al. 1997, Niazi et 
al. 2006, Razaq 2013), ectomycorrhizae have been described only for Russula 
brevipes Peck (Niazi et al. 2006). 

The aim of the present study is to elucidate the morphology and phylogeny of 
Russula sichuanensis and its ectomycorrhizae that form symbiotic associations 
with forest trees and play important role in the establishment and conservation 
of Himalayan forests. 

During a mushroom survey in 2012 and 2013, seven different collections of 
morphologically different basidiomata and their ectomycorrhizae were made 
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in Khyber Pakhtoonkhaw, Pakistan, which upon molecular analysis clustered 
together. Many of their morpho-anatomical features coincided with those of 
the Chinese species, Russula sichuanensis, and this identification was confirmed 
by ITS sequence analyses. Russula sichuanensis, a morphologically variable 
species, represents a new record for Pakistan, and its association with Pinus 
wallichiana has not previously been noted. Its ectomycorrhizae are described 
here for the first time. 


Materials & methods 


Morphological examination 

Basidiomata were collected in 2012 and 2013, photographed, vouchered, dried, 
and characterized morphologically. All specimens are conserved in the Herbarium, 
Department of Botany, University of the Punjab, Lahore, Pakistan (LAH). Sections of 
specimens were mounted in 5% KOH for microscopic observations under a biological 
microscope (MX4300H, Meiji Techno Co., Ltd., Japan). Tissues were mounted in 
phloxine to increase contrast of the structures, and Melzer’s reagent was used to examine 
amyloid reactions. 

Measurements of anatomical features (basidiospores, basidia, cystidia, pileus hyphae, 
and stipe hyphae) were made from at least 25 observations using an ocular micrometer 
and 100x oil-immersion objective; they include x = arithmetic mean of spore length 
and spore width for all spores measured, Q = spore length divided by spore width. Line 
drawings were made with a camera lucida. Color designations are from Munsell (1975). 


DNA extraction, PCR amplification, DNA sequencing 

Genomic DNAs were extracted from small pileus pieces using a modified CTAB 
method (Bruns 1995). Amplification of internal transcribed spacers (ITS1 and ITS2) 
and the 5.8S region of the nuclear ribosomal RNA gene were targeted using the primer 
pairs ITS1F/ITS4 (White et al. 1990; Gardes & Bruns. 1993) using the Extract-N-Amp 
plant DNA extraction Kit (Sigma-Aldrich, St. Louis, MO, USA). PCR was carried out 
under the following cycling parameters: initial denaturation (94°C for 1 min), 35 cycles 
(94°C for 1 min, 53°C for 1 min, and 72°C for 1 min), and final extension (72°C for 8 
min). Amplified PCR products were sent for purification and bidirectional sequencing 
to Macrogen (Korea). 


Sequence alignment & phylogenetic analysis 

Sequences were manually edited and assembled using BioEdit (www.mbio.ncsu.edu/ 
bioedit/bioedit.html). Following Dentinger et al. (2011) for complete ITS sequences, 
all sequences were trimmed with the conserved motifs 5’-( GAT) CATTA- and —GACCT 
(CAAA )-3’, and the alignment portions between them were included in analysis. All 
positions containing gaps and missing data were eliminated. There were a total of 406 
positions in the final dataset. The analysis involved 37 nucleotide sequences. Albatrellus 
ovinus was used as outgroup based on results reported by Li et al. (2013). 

Seven sequences from Russula sichuanensis basidiocarps, one sequence from an 
ectomycorrhizal root tip, and other related sequences retrieved from the GenBank 
were included in the phylogenetic analysis. The sequences were aligned by Muscle 


Russula sichuanensis ectomycorrhiza (Pakistan) ... 631 


Russula sichuanensis JX391968 
Hi Russula sichuanensis_KM596859 
Russula sichuanensis KF002768 
Russula sichuanensis KF002760 
H Russula sichuanensis_KM596860 
g1 | Russula sichuanensis JX391969, ex holotype 
E Russula sichuanensis_KM596861 
@ Russula sichuanensis_KM596862 
E Russula sichuanensis_KM596866 
97| |Æ Russula sichuanensis_KM596863 
E Russula sichuanensis_ KM596864 
E Russula sichuanensis_KM596865 
19 Russula_cessans_AY061730 
g2 | ı Russula_laricina_AY061685 
90 L Russula_nauseosa_JF908642 
Russula_curtipes AY061668 
Russula_violacea_AY061725 
63 Russula_caerulea_AY061661 
Russula_olivacea_AY061699 
Russula_decolorans_AY061670 
Russula_xerampelina_AY061734 
Russula_queletii_ AY061711 
Russula_exalbicans_AF418622 


Sect. Laricinae 


ail Russula_fellea_AY061676 
Russula_emetica_AY061673 
95 Russula_adusta_AY061652 
Russula_ingwa_JX266624 
53 Russula_delica_AY061671 


Russula7_foetens_AY061677 
Russula_cyanoxantha_AY061669 

Russula_ilicis AY061682 

96 Russula_crustosa_EU598194 
Russula_mustelina_AY061693 

77 Russula_aeruginea_AF418612 
79 Russula_grisea_AY061679 
Albatrellus_ovinus_AY198202 


FIG. 1. Phylogenetic relationship of Russula sichuanensis and closely related taxa based on Maximum 
Likelihood method inferred from nrITS sequences. The bootstrap consensus tree inferred from 
1000 replicates is taken to represent the evolutionary history of the taxa analyzed. 
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(Beta) using the default setting in Molecular Evolutionary Genetics Analysis (MEGA) 
software (Tamura et al. 2011). A phylogenetic tree was constructed with the Maximum 
Likelihood (ML) algorithm using a Jukes & Cantor (1969) model of nrITS sequences 
and nearest-neighbor-interchange (NNI) as ML heuristic search method using MEGA5 
software (Tamura et al. 2011). The topology was assessed by 1000 bootstrap replicates. 
Corresponding bootstrap values >50% are cited on the tree figure. Percent Identities 
(PID) and DNA divergences were calculated by DNAStar. 

The Russula sichuanensis sequences generated for this study were submitted to 
GenBank and the accession numbers for these and other closely related taxa used in the 
phylogenetic analysis are cited in phylogenetic tree (Fic. 1). 


Results 

Sequencing of the PCR products from the nrITS region of seven Pakistani 
specimens and one ectomycorrhizal root tip yielded fragments of 560-770 base 
pairs. Initial blast analysis of nucleotide sequences revealed a 99% maximum 
identity of the Pakistani sequences (GenBank KM596859-596866) with the 
sequence from the Russula sichuanensis holotype (GenBank JX391969). For 
phylogenetic analysis, several closely related taxa were selected from GenBank 
(FIG. 1). 

After removing the ambiguous letters from the aligned datasheet, 1008 
characters were subjected to phylogenetic analysis of which 294 characters 
were conserved, 428 were variable, 270 were parsimony informative, and 143 
were singletons. 

Percentage similarity was calculated by MegAlign (DNA Star Inc.). The 
Pakistani collections showed 98.6% identity and 0% genetic divergence with 
sequence JX391969 (R. sichuanensis holotype); 98% identity and 0% genetic 
divergence with sequence JX391968 (R. sichuanensis paratype), and 99% 
identity and 0.3% divergence with sequences KF002760 and KF002768. 


Taxonomy 


Russula sichuanensis G.J. Li & H.A. Wen, Mycotaxon 124: 179 (2013). Figs 2-4 

Basidiomata small to medium sized, epigeous or semi-hypogeous, russuloid 
when mature. PILEUS 15-62 mm diam., when young broadly parabolic, broadly 
convex with incurved margin and moderately to deeply indented centre or plane 
with flattened umbo, flesh grey upon drying, when mature straight or flaring 
to uplifted, tuberculate striate; surface dull, viscid, the centre very dark red 
(2.5R1/4) to very deep red (2.5R1/8), gradually less intense towards margin to 
dark red (2.5R3/4), the margin light grayish red (2.5R6/4) to deep red (2.5R3/10) 
or light olive brown (2.5Y5/2). LAMELLAE adnate or sinuate, regular, crowded, 
pale greenish yellow (7.5Y9/4) or light orange yellow (10YR9/8), margins even, 
concolorous. LAMELLULAE absent. STIPE 32-75 mm long, 7-11 um wide above, 
6-16 um wide below, central, clavate or subcylindrical, tapering upwards, dry, 


Russula sichuanensis ectomycorrhiza (Pakistan) ... 633 


NAY 


yo 
WZ 


1 D 


Fic. 2. Russula sichuanensis basidiomes: A, MSM0019; B, MSM0020; C, MSM0021; D, MSM0022. 
Scale bars-= 10 mm. 


dull, longitudinally rugulose, stuffed at first, becoming hollow (cavernous) 
upon maturity, white with yellowish tinge. Annulus and volva absent. Odor 
and taste not recorded. 

BASIDIOSPORES 8-12 x 7-10 um [x = 9 x 7.7 um, Q = 1.0-1.3], globose or 
subglobose, covered with amyloid warts linked as small crests and forming a 
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Fia. 3. Russula sichuanensis basidiomes: A, MSM0023; B, MSM0024. Scale bars = 10 mm. 


partial network, rarely intermixed with isolated verrucae, warts often 0.5-1.5 um 
high, thin-walled, yellowish in KOH, suprahilar plage well defined, 1.5-3.5 um. 
Basip1A 28-40 x 11-17 um, clavate, bi- or tetra-spored, thin-walled, hyaline 
in KOH; sterigmata (2.5-)4.5-6 um. TRAMA cellular, cells inflated, thin walled, 
hyaline or yellowish in KOH 19-36 x 22-47 um. CHEILOCYSTIDIA absent. 
PLEUROCYSTIDIA 42-78 x 8-16 um, clavate or rarely mucronate with crystal 
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Fic. 4. Russula sichuanensis (MSM0019): A, Basidia; B, Pleurocystidia; C, Pileocystidia; 
D, Caulocystidia; E, Basidiospores. Scale bars: A = 20 um; B-D = 25 um; E = 8 um. 


or refringent contents, hyaline or yellowish in KOH, thin walled. PILEIPELLIS 
an ixotrichoderm composed of cylindrical hyphae, 3.5-6.3 um diam., thin- 
walled, hyaline or greenish in KOH. PILEOCYSTIDIA with refractive contents, 
septate, clavate to cylindrical, 6.3-11.3 um. STIPITIPELLIS interwoven, hyphae 
cylindrical, 3.3-5.3 um, septate, hyaline in KOH. CauLocystip1a with 
refractive contents, septate, clavate to cylindrical, 56-89 x 6-10 um. CLAMP 
CONNECTIONS absent. 
MATERIAL EXAMINED (all in LAH): PAKISTAN, KHYBER PAKHTOON KHaAw, 
Mansehra, Chattar plain, under Pinus wallichiana A.B. Jacks., 15 September 2012, 
M. Saba MSM0019 (GenBank KM596859), MSM0020 (GenBank KM596860); 22 
September 2013, M. Saba, MSM0021 (GenBank KM596861), MSM0022 (GenBank 
KM596862), MSM0023 (GenBank KM596863), MSM0024 (GenBank KM596864), 
MSM0025 (GenBank KM596865); Shangla, ectomycorrhiza near rhizosphere of Pinus 
wallichiana, 2 September 2013, M. Saba MSE003 (GenBank KM596866). 


Morphology of the ectomycorrhizae Fic. 5 

ECTOMYCORRHIZAL SYSTEM: Dichotomous, ramified system 2.5 mm long. 
Unramified ends straight, 0.5 mm long, orange-brown to grayish brown and 
0.5 mm in diameter (basal portion 0.4 mm wide; apical portion 0.5 mm wide). 
Very tips orange brown. Mantle surfaces smooth, matte, rarely giving rise to 
emanating hyphae and with host visible through mantle. Mycorrhizae orange- 
brown when young and grayish brown when old. 
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Fic. 5. Russula sichuanensis (MSE003) ectomycorrhizae: A, Morphotypes; B, Outer mantle; 
C, Inner mantle; D, Emanating hyphae. Scale bars: A = 5 mm; B-D = 15 um. 


OUTER MANTLE: pseudoparenchymatous (type L; Agerer 1987-2012) with 
irregularly shaped, often interlocking cells. Cells olive yellow to yellowish 
brown, 9.6-28.3 x 6.3-17.4 um. INNER MANTLE: pseudoparenchymatous (type 
L; Agerer 1987-2012) with irregularly shaped, often interlocking cells. Cells 
olive yellow to yellowish brown, 47-86.5 x 24-58 um. EMANATING HYPHAE: 
curved, frequently septate, yellowish brown to olive yellow, thin walled, rarely 
branched, 8-25.9 x 2.4-4.1 um. No clamp connection observed; H-type 
anastomosis observed among hyphae. RHIZOMORPHS: not observed. CYSTIDIA: 
not observed. 
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Discussion 

Russula sichuanensis was originally described by Li et al. (2013) from Sichuan 
province, China. It is characterized by an uneven pileus that often covers a 
white percurrent stipe and tightly crowded sinuate cream to yellow lamellae; 
it was additionally described as having a cream-tinged pileus (often olivaceous 
brown at the centre), a yellowish to orange sublamellate gleba, and stipitate 
basidiomata (Li et al. 2013). Our collections differ morphologically from the 
type description as follows: (1) curly and tightly contorted pilei were observed 
in only a few young basidiomata while the majority of mature basidiomata 
have broadly convex pilei with straight or uplifted margins and regular, straight 
(without convolutions) lamellae; (2) the pileus color is variable, ranging from a 
very dark red to very deep red centre, dark red gradually towards margin, light 
grayish red to deep red margin, or light olive brown. Most of the microscopic 
features, including shape, size of basidiospores, their ornamentation, 
presence of ixotrichoderm pileipellis with septate pileocystidia and septate 
caulocystidia, coincide with the description given by Li et al. (2013); however 
our collections had basidiospores with a well defined suprahilar plage and 
frequent pleurocystidia. Nonetheless, a BLAST comparison of all ITS sequences 
retrieved from our seven different basidiomata and one ectomycorrhizal root 
tip with the holotype of R. sichuanensis showed 99% identity and 100% query 
coverage and clustered with the holotype sequence and three other Chinese 
sequences in a clade with 91% support in phylogenetic reconstruction (Fic. 1). 

Russula sichuanensis can now be characterized as a morphologically variable 
species, especially with respect to pileus colour and shape and which produces 
either a sublamellate gleba or regular lamellae. In addition, it can be characterized 
by partially small to medium sized stipitate basidiomata, subreticulate globose 
to subglobose basidiospores, and pluriseptate pileocystidia. The morphological 
variation in colour and forms result from non-genetic phenotypic plasticity 
and some edaphic or substrate factors may be involved in maintaining natural 
variants (Nicholl 1996; Petersen & Hughes 1999). Phenotypic plasticity in 
mushrooms has been largely unexplored, no doubt due to the ephemeral nature 
of most fungal basidiomata, which form only when environmental conditions 
are appropriate and then last only a short time. 

Our phylogeny clusters R. sichuanensis with members of Russula sect. 
Laricinae, suggesting that it is closely related to the R. laricina/cessans group 
as previously depicted by Li et al. (2013). Typical morphological characters 
of this section include a cavernous stipe, tuberculate-striate cap margin, and 
flesh that does not yellow but greys after drying, all characters observed in our 
collections. 

Russula mycorrhizae can be characterized by a pseudoparenchymatous 
outer mantle with either angular or irregularly shaped cells and which may be 
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more or less densely covered with cystidia or acystidiate (Eberhardt 2002). In 
addition to the mantle features, Agerer (1999) defined a russuloid rhizomorph 
type, distinguished by its appression to a substrate (often to thicker roots) 
and differentiated into a peripheral layer with or without cystidia and central 
layer with vessel-like hyphae. According to Eberhardt (2002), ectomycorrhizae 
in R. subg. Russula, Incrustatula, Tenellula (in which R. sichuanensis falls), 
Polychromidia, and Insidiosula and in R. sect. Felleinae are acystidiate with a 
pseudoparenchymatous outer mantle that has irregularly shaped and often 
interlocking cells. 

The ectomycorrhizal characters of R. sichuanensis are same as those of Russula 
subg. Tenellula, in which Li et al. (2013) described this species: dichotomously 
branched and orangish brown to grayish brown tips, a pseudoparenchymatous 
and acystidiate mantle with irregularly shaped and often interlocking cells and 
frequently septate, curved emanating hyphae. The presence of rhizomorphs 
was not noted. 

Li et al. (2013) reported R. sichuanensis from conifer forests (Picea spp.); 
our study describes the basidiocarp and its ectomycorrhiza from monoculture 
Pinus wallichiana forests located in Himalayan moist temperate region. 
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